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Abstract. The ATLAS liquid argon calorimeter has been operating continuously since August 2006. At 
this time, only part of the calorimeter was readout, but since the beginning of 2008, all calorimeter cells 
have been connected to the ATLAS readout system in preparation for LHC collisions. This paper gives an 
overview of the liquid argon calorimeter performance measured in situ with random triggers, calibration 
data, cosmic muons, and LHC beam splash events. Results on the detector operation, timing performance, 
electronics noise, and gain stability are presented. High energy deposits from radiative cosmic nmons and 
beam splash events allow to check the intrinsic constant term of the energy resolution. The uniformity of 
the electromagnetic barrel calorimeter response along r] (averaged over (j)) is measured at the percent level 
using minimum ionizing cosmic muons. Finally, studies of electromagnetic showers from radiative muons 
have been used to cross-check the Monte Carlo simulation. The performance results obtained using the 
ATLAS readout, data acquisition, and reconstruction software indicate that the liquid argon calorimeter 
is well-prepared for collisions at the dawn of the LHC era. 
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1 Introduction 

Installation of the liquid argon (LAr) calorimeter in the 
ATLAS [T] experimental hall was completed in early 2008. 
Until recently, the expected performance of the LAr cal- 
orimeter was extrapolated from intensive testing of a few 
modules with electron and pion beams from 1998 to 2003 
(Ref. [ H ll l li l lS l lg llTllSll^fTU] ] . and in 2004 of a complete AT- 
LAS detector slice PTIIT^ITS] . The 20 months separating 
the completion of the installation from the first LHC col- 
lisions have been used to commission the LAr calorim- 
eter. This paper reviews the first in situ measurements 
of the electronics stability, the quality of the energy re- 



construction, the calorimeter response uniformity and the 
agreement between data and the Monte Carlo simulation 
of electromagnetic shower shapes. The measurements are 
performed using calibration triggers, cosmic muons, and 
the first LHC beam events collected during this 20 months 
period. The results and the experience gained in the oper- 
ation of the LAr calorimeter provide the foundation for a 
more rapid understanding of the experimental signatures 
of the first LHC collisions, involving electrons, photons, 
missing transverse energy (E!^^^^), jets, and rs where the 
LAr calorimeter plays a central role. 
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This paper is organized as follows. Section [2] gives the 
present hardware status of the LAr calorimeter. Section [3] 
details the level of understanding of the ingredients enter- 
ing the cell energy reconstruction: pedestals, noise, elec- 
tronic gains, timing, and the quality of the signal pulse 
shape predictions. The current understanding of the first 
level trigger energy computation is also discussed. Sec- 
tion 2] describes the in situ performance of the electromag- 
netic LAr calorimeter using ionizing and radiating cosmic 
muons. Lastly, Section [5] draws the conclusions. 



2 LAr calorimeter hardware status and data 
taking conditions 




LAr forward (FCal) 



The LAr calorimeter is composed of electromagnetic and 
hadronic sub-detectors of which the main characteristics 
are described in Section [^H During the detector and elec- 
tronics construction and installation, regular and stringent 
quality tests were performed, resulting in a fully functional 
LAr calorimeter. The operational stability of the cryostats 
since March 2008 is discussed in Section [221 The current 
status of the high voltage and the cell readout are dis- 
cussed in Sections 12.31 and 12.41 respectively. Finally, the 
general data taking conditions are given in Section [2.5l In 
ATLAS, the positive x-axis is defined as pointing from the 
interaction point to the center of the LHC ring, the pos- 
itive y-axis is defined as pointing upwards, and the posi- 
tive z-axis corresponds to protons running anti-clockwise. 
The polar angle 9 is measured from the beam axis (z- 
axis) , the azimuthal angle (j) is measured in the transverse 
(a;?/)-plane, and the pseudorapidity is defined aa rj = -In 
tan(6l/2). 



2.1 Main characteristics of the LAr calorimeter 

The LAr calorimeter [T], shown in Figure [TJ is composed of 
sampling detectors with full azimuthal symmetry, housed 
in one barrel and two endcap cryostats. More specifically, 
a highly granular electromagnetic (EM) calorimeter with 
accordion-shaped electrodes and lead absorbers in liquid 
argon covers the pseudorapidity range |?7| < 3.2, and con- 
tains a barrel part (EMB [T3], |?7| < 1.475) and an endcap 
part (EMEC [H], 1.375 < |7y| < 3.2). For \t]\ < 1.8, a pre- 
sampler (PS [TBlfTS] ). consisting of an active LAr layer 
and installed directly in front of the EM calorimeters, 
provides a measurement of the energy lost upstream. Lo- 
cated behind the EMEC is a copper-liquid argon hadronic 
endcap calorimeter (HEC [17], 1.5 < \r]\ < 3.2), and a 
copper/tungsten-liquid argon forward calorimeter (FCal 
[TB]) covers the region closest to the beam at 3.1 < \r]\ < 
4.9. An hadronic Tile calorimeter (|r;| < 1.7) surrounding 
the LAr cryostats completes the ATLAS calorimetry. 

All the LAr detectors are segmented transversally and 
divided in three or four layers in depth, and correspond 
to a total of 182,468 readout cells, i.e. 97.2% of the full 
ATLAS calorimeter readout. 



Fig. 1. Cut-away view of the LAr calorimeter, 17 m long (bar- 
rel -f endcaps) and 4 rn of diameter. 



The relative energy resolution of the LAr calorimeter 
is usually parameterized by: 



E 



E 



(1) 



where (a) is the stochastic term, (6) the noise term and 
(c) the constant term. The target values for these terms 
are respectively a ~ 10%, h ~ 170 MeV (without pile-up) 
and c = 0.7%. 



2.2 Cryostat operation 

Variations of the liquid argon temperature have a direct 
impact on the readout signal, and consequently on the en- 
ergy scale, partly through the effect on the argon density, 
but mostly through the effect on the ionization electron 
drift velocity in the LAr. Overall, a — 2%/K signal varia- 
tion is expected [TH]. The need to keep the corresponding 
contribution to the constant term of the energy resolution 
(Eq. [1]) negligible (i.e. well below 0.2%) imposes a tem- 
perature uniformity requirement of better than 100 mK 
in each cryostat. In the liquid, ^-^500 temperature probes 
(PTIOO platinum resistors) are fixed on the LAr detec- 
tor components and read out every minute. In 2008-2009, 
installation activities in the ATLAS cavern prevented a 
stable cryostat temperature. A quiet period of ten days 
around the 2008 Christmas break, representative of what 
is expected during LHC collisions, allowed a check of the 
temperature stability in the absence of these external fac- 
tors. The average dispersion (RMS) of the measurements 
of each temperature probe over this period is 1.6 mK (5 
mK maximum) , showing that no significant local temper- 
ature variation in time is observed in the three cryostats. 
Over this period, the temperature uniformity (RMS of all 
probes per cyostat) is illustrated for the barrel in Figured] 
and gives 59 mK. Results for the two endcap cryostats are 
also in the range 50-70 mK, below the required level of 
100 mK. The average cryostat temperatures are slightly 
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different for the barrel (88.49 K) and the two endcaps 
(88.67 and 88.45 K) because they are independently reg- 
ulated. An energy scale correction per cryostat will there- 
fore be applied. 



EM BARREL 




7 88.8 88.9 
Temperature [K] 

Fig. 2. Distribution of barrel cryostat probe temperatures av- 
eraged over a period of ten days. 



To measure the effects of possible out-gassing of calo- 
rimeter materials under irradiation, which has been mini- 
mized by careful screening of components, 30 purity moni- 
tors measuring the energy deposition of radioactive sources 
in the LAr are installed in each cryostat and read every 
15 minutes. The contribution to the constant term of the 
energy resolution is negligible for a level of electronega- 
tive impurities below 1000 ppb O2 equivalent. All argon 
purity measurements over a period of two years are stable, 
in the range 200 ± 100 ppb O2 equivalent, well below this 
requirement. 

In summary, measurements of the liquid argon tem- 
perature and purity demonstrate that the stability of the 
operation of the three LAr cryostats is in the absence of 
proton beams within the required limits ensuring a negli- 
gible contribution to the energy resolution constant term. 



region. For example, in the EM calorimeter, faulty elec- 
trodes were connected to separate HV lines during the 
assembly phase at room temperature while, if the defect 
was identified during cryostat cold testing, the high volt- 
age sector was divided into two in (j), each connected sepa- 
rately. The effect of zero voltage on one side of an electrode 
was studied in beam tests proving that with offline correc- 
tions the energy can still be measured, with only a small 
loss in accuracy. Finally, for HV sectors with a permanent 
short-circuit, high voltage modules permitting large DC 
current draws of up to 3 mA (more than three orders of 
magnitude above the nominal limit) are used in order to 
operate the faulty sector at 1000 V or above. 

As a result, 93.9% of readout cells are operating under 
nominal conditions and the rest sees a reduced high volt- 
age. However, even with a reduced high voltage, signals 
can be well reconstructed by using a correction scale fac- 
tor. Figure [3] shows the distribution of all HV correction 
factors for the EM, HEC and FCal cells as of the end of 
September 2009. Since the beginning of 2008, no changes 
have been observed. The largest correction occurs if one 
side of an EM electrode is not powered, and only half of 
the signal is collected. For the faulty cells, this correction 
factor is applied online at the energy reconstruction level. 
A similar correction is currently being implemented at the 
first level (LI) trigger. 




1.6 1.8 2 
HV correction factor 



2.3 High voltage status 

The electron/ion drift speed in the LAr gap depends on 
the electric field, typically 1 kV/mm. Sub-detector-specific 
high voltage (HV) settings are applied. In the EM barrel, 
the high voltage is constant along 77, while in the EMEC, 
where the gap varies continuously with radius, it is ad- 
justed in steps along rj. The HV supply granularity is typ- 
ically in sectors of Ar] x A(j) = 0.2 x 0.2. For redundancy, 
each side of an EM electrode, which is in the middle of the 
LAr gap, is powered separately. In the HEC, each sub-gap 
is serviced by one of four different HV lines, while for the 
FCal each of the four electrode groups forming a normal 
readout channel is served by an independent HV line. 

For HV sectors with non-optimal behavior, solutions 
were implemented in order to recover the corresponding 



Fig. 3. High voltage correction factors for all LAr cells at the 
end of September 2009. 

In conclusion, since the beginning of 2008, all 182,468 
readout cells are powered with high voltage, and no dead 
region exists. Signals from regions with non-nominal high 
voltage are easily corrected and their impact on physics is 
negligible. 

2.4 Readout cell status 

The cell signals are read out through 1524 Front-End 
Boards (FEBs [201121] ') with 128 channels each, which sit 
inside front-end crates that are located around the periph- 
ery of the cryostats. The FEBs perform analog processing 



15 



(amplification and shaping - except for tlie HEC wliere tlie 
amplification is done inside the cryostat), store the signal 
while waiting for the LI trigger decision, and digitize the 
accepted signals. The FEBs also perform fast analog sum- 
ming of cell signals in predefined projective "towers" for 
the LI trigger. 

The digitized signals are transmitted via optical fibers 
to the Readout Drivers (RODs) t22j located in the count- 
ing room 70 m away. The cell energy is reconstructed on- 
line in the ROD modules up to a nominal maximum LI 
rate of 75 kHz. The cell and trigger tower energy recon- 
struction is described in detail in Section [31 

The response of the 182,468 readout cells is regularly 
monitored using 122 calibration boards (33] located in the 
front-end crates. These boards inject calibrated current 
pulses through high-precision resistors to simulate energy 
deposits in the calorimeters. At the end of September 
2009, 1.3% of cells have problems. The majority of them, 
i.e. 1.2% of the total number of cells, are not read-out 
because they are connected to 17 non-functioning FEBs. 
On these FEBs, the active part (VCSEL) of the optical 
transmitter to the ROD has failed. This failure, occur- 
ring at a rate of two or three devices per month, is un- 
der intensive investigation and are expected to be fixed 
during the next LHC shutdown. The remaining 0.1% of 
cells with problems can be split in three sub- types: in- 
curable cells, i.e. cells not responding to the input pulse 
(0.02%), or which are permanently (0.03%) or sporadi- 
cally (0.07%) very noisy. The first two types are always 
masked in the event reconstruction (121 cells), while the 
sporadically very noisy cells, not yet well understood, are 
masked on an event by event basis. For cells which do 
not receive calibration signals (0.3%) average calibration 
constants computed among neighboring cells are used. For 
cells with non-nominal high voltage (6.1%) a software cor- 
rection factor is applied. Both have very limited impact on 
the energy reconstruction. 

In total, 180,128 cells, representing 98.7% of the total 
number of cells in the LAr calorimeter, are used for event 
reconstruction at the end of September 2009. The number 
of inactive cells (1.3%) is dominated by the cells lost due to 
faulty optical drivers (1.2%): apart from these, the number 
of inactive cells has been stable in time. 



2.5 Data taking conditions 

The results presented here focus on the period starting in 
September 2008 when all the ATLAS sub-detectors were 
completed and integrated into the data acquisition. Apart 
from regular electronics calibration runs, two interesting 
types of data are used to commission the LAr calorime- 
ter: the beam splash events and the cosmic muons. The 
first type corresponds to LHC events of September lO"^ 
2008 when the first LHC beam hit the collimators located 
200 m upstream of the ATLAS interaction point. A cas- 
cade of pions and muons parallel to the beam axis fired 
the beam related trigger, illuminated the whole ATLAS 
detector and deposited several PeV per event in the LAr 
calorimeter. The second type corresponds to long cosmic 



muon runs acquired on September-October 2008 and on 
June- July 2009 where more than 300 million events were 
recorded, corresponding to more than 500 TB of data. 

For the LAr commissioning, LI calorimeter triggers are 
used to record radiative energy losses from cosmic muons 
while the first level muon spectrometer and second level 
inner detector triggers are used to study pseudo-pro jective 
minimum ionizing muons. In most of the runs analyzed, 
the toroidal and solenoidal magnetic fields were at the 
nominal value. 



3 Electronic performance and quality of cell 
energy reconstruction 

The robustness of the LAr calorimeter energy reconstruc- 
tion has been studied in detail using calibration and ran- 
domly triggered events, cosmic muons and beam splash 
events. Section lHTTI briefiv describes the energy reconstruc- 
tion method in the trigger towers and in the cells, as well 
as a validation study of the trigger. The time stability of 
the electronics is discussed in Section[X2J The status of the 
electronics timing for the first LHC collisions is presented 
in Section 13.31 and the quality of the LAr calorimeter en- 
ergy reconstruction is assessed in Section 13.41 

3.1 Energy reconstruction in the LAr calorimeter 

When charged particles cross the LAr gap between elec- 
trodes and absorbers, they ionize the liquid argon. Under 
the infiuence of the electric field, the ionization electrons 
drift towards the electrode inducing a current. The ini- 
tial current is proportional to the energy deposited in the 
liquid argon. The calorimeter signals are then used to com- 
pute the energy per trigger tower or per cell as discussed 
in this section. 



3.1.1 Energy reconstruction at the first level calorimeter 
trigger 

The timing requirements for the LI trigger latency can 
only be met with fast analogue summing in coarse gran- 
ularity. In the EM part, the pre-summation of analog 
signals per layer on the FEBs serves as input to tower 
builder boards where the final trigger tower signal sum 
and shaping is performed. In the HEC and FCal, the 
summation is performed on the FEBs and transmitted to 
tower driver boards where only shaping is done. The tower 
sizes are Ar] x Acf) = 0.1 x 0.1 for \ri\ < 2.5 and go up to 
Aijx Acf)^ 0.4 X 0.4 for 3.1 < \ri\ < 4.9. The analog trig- 
ger sum signals are sent to receiver modules in the service 
cavern. The main function of these modules is to compen- 
sate for the differences in energy calibration and signal 
attenuation over the long cables using programmable am- 
plifier gains ((7r). The outputs are sent to LI trigger pre- 
processor boards which perform the sampling at 40 MHz 
and the digitization of five samples. At this stage, both 
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the transverse energy and bunch crossing are determined 
using a finite impulse response filter, in order to maximize 
the signal-to-noise ratio and bunch crossing identification 
efficiency. During ATLAS operation, the output QrA^^ of 
the filter, which uses optimal filtering, is passed to a look- 
up table where pedestal (P in ADC counts) subtraction, 
noise suppression and the conversion from ADC counts to 
transverse energy in GeV (pADC-s-Gev) performed in or- 
der to extract the final transverse energy value (EJ^^) for 
each trigger tower: 

Ek' = ^^ADC^Gev(5R^" - 9rP)- (2) 

Arrays (in rj—cf)) of these EJ^^ energies, merged with similar 
information coming from the Tile calorimeter, are subse- 
quently used to trigger on electrons, photons, jets, rs and 
events with large missing transverse energy. 

3.1.2 Energy reconstruction at cell level 

At the cell level, the treatment of the analog signal is also 
performed in the front-end electronics. After shaping, the 
signal is sampled at 40 MHz and digitized if the event 
was selected by the LI trigger. The reconstruction of the 
cell energy, performed in the ROD, is based on an opti- 
mal filtering algorithm applied to the samples Sj j24j . The 
amplitude A, in ADC counts, is computed as: 

i=i 

where p is the ADC pedestal fSection I3.2.ip . The Opti- 
mal Filtering Coefficients (OFCs) aj are computed per 
cell from the predicted ionization pulse shape and the 
measured noise autocorrelation to minimize the noise and 
pile-up contributions to A. For cells with sufficient signal, 
the difference [At in ns) between the digitization time and 
the chosen phase is obtained from: 

^ (4) 

i=i 

where hj are time-OFCs. For a perfectly timed detector 
and in-time particles \At\ must be close to zero, while 
larger values indicate the need for better timing or the 
presence of out-of-time particles in the event. 

The default number of samples used for A and At 
computation is Ngampios = 5, but for some specific analy- 
ses more samples, up to a maximum of 32, are recorded. 
Finally, including the relevant electronic calibration con- 
stants, the deposited energy (in MeV) is extracted with: 

£^ccll = PpA^MoV X PdAC^mA X Mphys X G X A, (5) 

Mcali 

where the various constants are linked to the calibration 
system: the cell gain G (to cover energies ranging from a 
maximum of 3 TeV down to noise level, three linear gains 



are used: low, medium and high with ratios ^ 1/10/100) 
is computed by injecting a known calibration signal and 
reconstructing the corresponding cell response; the fac- 
tor 1 / M c aii'^ quantifies the ratio of response to a calibra- 
tion pulse and an ionization pulse corresponding to the 
same input current; the factor Pdac->-aiA converts digital- 
to-analog converter (DAC) counts set on the calibration 
board to ^A] finally, the factor P^A^-MeV is estimated 
from simulations and beam test results, and includes high 
voltage corrections for non- nominal settings (see Sec l2.3p . 
Note that the crosstalk bias in the finely segmented first 
layer of the electromagnetic calorimeter is corrected for in 
the gain G g]. 

3.1.3 Check of the first level tower trigger energy 
computation 

The trigger decision is of utmost importance for ATLAS 
during LHC collisions since the data-taking rate is at max- 
imum 200 Hz because of bandwidth limitations, i.e. a fac- 
tor 2 X 10^ smaller than the 40 MHz LHC clock. It is 
therefore important to check that no systematic bias is 
introduced in the computation of the LI trigger energy 
and that the trigger energy resolution is not too degraded 
with respect to the offline reconstruction. In the following, 
this check is performed with the most granular part of the 
LAr calorimeter, the barrel part of the EM calorimeter, 
where 60 cell signals are summed per trigger tower. 

Since cosmic muon events occur asynchronously with 
respect to the LHC clock, and the electronics for both the 
trigger and the standard readout is loaded with one set of 
filtering coefficients (corresponding to beam crossing), the 
reconstructed energy is biased by up to 10%, depending 
on the phase. For the study presented here, A^^ is recom- 
puted offline by fitting a second-order polynomial to the 
three highest samples transmitted through the processors. 
The most critical part in the trigger energy computation is 
then to calibrate the individual receiver gains g^i- For that 
purpose, a common linearly increasing calibration pulse is 
sent to both the LI trigger and the normal cell circuits: 
the inverse receiver gain 1/(?r is obtained by fitting the 
correlation between the LI calorimeter transverse energy 
{E\^) and the sum of cell transverse energies in the same 
trigger tower, later called ojfline trigger tower {E}^^^). In 
cosmic muon runs, receiver gains are set to 1.0 and are 
recomputed offline with dedicated calibration runs. As a 
cross check, the gain was also extracted using LHC beam 
splash event data which covers the full detector. In both 
cases, the LI transverse energy is computed as in Eq. [2j 

In the EM calorimeter, radiating cosmic muons may 
produce a local energy deposit of a few GeV, and fire 
the EM calorimeter trigger condition EM3 that requires 
a transverse energy greater than 3 GeV in a sum of four 
adjacent EM trigger towers. To mimic an electron coming 
from the interaction point, only those events that contain 
a track reconstructed with strict projectivity cuts are con- 
sidered. Here, the LI calorimeter transverse energy is com- 
puted using the gains determined with calibration runs. 
Figure S] shows the correlation between E^^ and E^^'^ . 
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Computing the ratio of E^^ and E^^'^ gives a Gaussian 
distribution with a mean of 1.015±0.002, showing the very 
good correspondence between these two quantities, espe- 
cially at low energy. This also shows that the trigger en- 
ergy is well calibrated and almost unbiased with respect 
to the LAr readout. 
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Fig. 4. LI transverse energy (E^^ ) computed with the receiver 
gains extracted from calibration runs versus the sum of cell 
transverse energies in the same trigger tower (E^^'' ). 



As a crosscheck, a similar study was performed with 
gains computed from the beam splash events, without the 
projectivity cut. A slight degradation of the resolution is 
observed at high energy, but not at low energy where the 
noise dominates. Taking advantage of the higher statistics, 
it is possible to compute the 5 GeV "turn-on curve" , i.e. 
the relative efficiency for an offline trigger tower to meet 
the requirement E^^ > 5 GeV as a function of E^^'' . This 
is not the absolute efficiency as the calorimeter trigger 
condition EMS is used to trigger the events. The efficiency 
is shown in Figure [Bl where a sharp variation around a 
EI^^ = 5 GeV energy threshold is observed. 
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Figure [5] shows the corresponding resolution computed 
as the relative difference of i?^^ and EJ^^^ . At low energy, 
the difference is dominated by electronic noise since the 
two readout paths have only part of their electronics in 
common. The ATLAS specification of 5% of LI transverse 
energy resolution is reached for energies greater than 10 
GeV. The LI transverse energy resolution reaches around 
3% at high energy. 



Fig. 6. Turn-on curve efficiency for E^^ > ^ GeV requirement 
obtained with events triggered by the EM3 LI Calorimeter trig- 
ger. 

These results give confidence that EM showers (elec- 
trons and photons) will be triggered efficiently in LHC 
events. After this study, the gains gn were extracted from 
dedicated calibration runs and loaded into the receivers to 
be used for the first LHC collisions. 
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3.2 Electronic stability 

Hundreds of millions of randomly triggered and calibra- 
tion events can be used for a study of the stability of the 
properties of each readout channel, such as the pedestal, 
noise and gain. The first two quantities are computed for 
each cell as the mean (pedestal) of the signal samples Sj 
in ADC counts, and the width (noise) of the energy distri- 
bution. The gain is extracted by fitting the output pulse 
amplitudes against calibration pulses with increasing am- 
plitudes. 



3.2.1 Pedestal 



Fig. 5. Relative difference of E^^ and E^^'^ (LI Calorimeter 
Et resolution) as a function of E^^^ . Strict projectivity cuts 
for the track pointing to the EM shower are applied. Horizontal 
error bars reflect the RMS of E^^'' in each bin. 



The stability of the pedestals is monitored by measuring 
variations with respect to a reference pedestal value for 
each cell. For each FEB, an average over the 128 channels 
is computed. 
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As an example, Figure [7] illustrates the results for the 
48 HEC FEBs over a period of six months in 2009. A 
slight drift of the pedestal with time, uncorrelated with the 
FEB temperature and/or magnetic field configurations, is 
observed. Overall, the FEB pedestal variations follow a 
Gaussian distribution with a standard deviation of 0.02 
ADC counts, i.e. below 2 MeV. The same checks have been 
performed on all other FEBs, and give typical variations of 
around 1 (0.1) MeV and 10 (1) MeV in the EM and FCal 
calorimeters respectively, in medium (high) gain. These 
variations are much lower for the EM and HEC or at the 
same level for the FCal than the numerical precision of 
the energy computation, which is 8 (1) MeV in medium 
(high) gain. 
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Fig. 7. Average FEB pedestal variations in ADC counts, in 
medium gam, for the HEC during 6 months of data taking in 
2009. The crosses indicate the mean value for each ttme sltce. 



Fig. 8. Electronic noise (cTnoiso^ in randomly triggered events at 
the EM scale in individual cells for each layer of the calorimeter 
as a function of \r]\. Results are averaged over 4>. 

coherent noise has been estimated to 2%, by studying si- 
multaneous increase of noise in a group of channels. 

Systematic studies of noise stability have been pur- 
sued: all noise variations are typically within ± 1 keV, 0.1 
MeV and 1 MeV for EM, HEC and FCal, respectively. No 
correlations with the FEB temperature and/or changes of 
magnetic field conditions have been observed. 



3.2.3 Gain 



During the LHC running, it is foreseen to acquire pedesi 
and calibration runs between fills, thus it will be possible 
to correct for any small time dependence such as observed 
in Figure [71 In the same spirit, random triggers collected 
during physics runs can be used to track any pedestal 
variations during an LHC fill. 

3.2.2 Noise 

Figure |5] shows the noise measured in randomly triggered 
events at the cell level as a function of 77 for all layers 
of the LAr calorimeters. In all layers, a good agreement 
with the expected noise is observed. Noise values are 
symmetric with respect to 77 = and uniformly in (p within 
few percents. In the EM calorimeters, the noise ranges 
from 10 to 50 MeV, while it is typically a factor of 10 
greater in the hadronic endcap and forward calorimeters 
where the granularity is 20 times coarser and the sampling 
fractions are lower. It should be noted that these results 
are obtained using five samples in Eq. ^ and ([5]), i.e. 
the noise is reduced by a factor varying from 1.5 to 1.8, 
depending on -q, with respect to the single-sample noise 
value. 

The coherent noise over the many cells used to mea- 
sure electron and photon energies in the EM calorimeters 
should be kept below 5% [25 of the incoherent noise (i.e. 
the quadratic sum of all channel noise). For the second 
layer of the EM calorimeter, the contribution from the 



The calibration pulse is an exponential signal (controlled 
by two parameters, /stop and Tcaii) which emulates the 
triangular ionization signal. It is injected on the detec- 
tor as close as possible to the electrodes, except for the 
FCal where it is applied at the base-plane of the front-end 
crates [TH]- Thus, the analog cell response is treated by 
the FEBs in the same way as an ionization signal, but it 
is typically averaged over 100 triggers in the RODs and 
transmitted offline where the average signal peak height 
is computed. The cell gain is extracted as the inverse ra- 
tio of the response signal in ADC counts to the injected 
calibration signal in DAC counts. 

The stability of the cell gain is monitored by looking at 
the relative gain difference averaged over 128 FEB chan- 
nels. This is illustrated in Figure |9] for the 1448 FEBs of 
the EM calorimeter, in high gain. All variations are within 
±0.3% and similar results are obtained for medium and 
low gains. An effect of 0.2% on the gains has recently been 
identified as coming from a particular setting of the FEBs. 
The two populations are most probably coming from this 
effect. Regular update of calibration database take ac- 
count of the variations. Similar results are obtained for 
the HEC, and variations within ±0.1% are measured for 
the FCal. 

In conclusion, results presented for the pedestals, noise, 
and gains illustrate the stability of the LAr electronics 
over several months of data taking. Values are stored in 
the ATLAS calibration database and are used for online 
and offline reconstruction. 
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Fig. 9. Average FEB (high) gain variations during 6 months 
of 2009 data taking, m the EM part of the calorimeter. The 
crosses indicate the mean value for each time slice. 

3.2.4 Global check with E^^'"''' variable 

Another way to investigate the level of understanding of 
pedestals and noise in the LAr calorimeter is to compute 
global quantities in randomly triggered events with the 
calorimeter, such as the vector sum of transverse cell en- 
ergies. The calorimetric missing transverse energy E^^^^ is 
defined as: 

_gmiss ^ _^;^-"S^sin0,COS0„ 

i^r" = -Y^^IE^^^^^^^,, (6) 

where Ei is the cell energy, 9i its polar angle and 4>i its az- 
imuthal angle. Because of the high granularity of the LAr 
calorimeter, it is crucial to suppress noise contributions 
to E!^^^^, i.e. limit the number of cells, iVceii, used in the 
sum. In ATLAS, this is done with two methods: i) a cell- 
based method in which only cells above a noise threshold 
of two standard deviations {\Ei\ > 2crnoisc) are kept; ii) a 
cluster-based method which uses only cells belonging to 
three-dimensional topological clusters [55] . These clusters 
are built around \Ei\ > 4CTnoiso seeds by iteratively gath- 
ering neighboring cells with \Ei\ > 2(Tnoiso and, in a final 
step, adding all direct neighbors of these accumulated sec- 
ondary cells (Topocluster 4/2/0). In randomly triggered 
events, about 8500 and 500 LAr cells, respectively, are 
selected with these two noise-suppression methods. 

The distributions of and E"'"'''' should be Gaus- 

sian and centered on zero in randomly triggered events. 
The measurements are compared with a Gaussian noise 
model, where no pedestal shift or coherent noise is present, 
obtained by randomizing the cell energy according to a 
Gaussian model for the cell noise. For this E^^^^ com- 
putation, cells with very high noise (see Section 12. 4|) are 
removed from the computation. 

Figure [TO] shows the E!^^^^ distributions for a randomly 
triggered data sample acquired in 15 hours. The two noise 
suppression methods are compared to the corresponding 
Gaussian noise model. For the cell-based method, a good 
agreement is observed between the data and the simple 
model. Because of the lower number of cells kept in the 
cluster-based method, a smaller noise contribution to E!^^^^ 
is observed. The agreement between the data and the 



model is not as good as for the cell-based method, re- 
flecting the higher sensitivity of the cluster-based method 
to the noise description. In both cases, no E^^^^ tails are 
present, reflecting the absence of large systematic pedestal 
shifts or abnormal noise. 



RANDOM TRIGGER LAr CALORIMETER 




Fig. 10. E^^^^ distribution with LAr calorimeter cells for 
135,000 randomly triggered events in June 2009. The dots 
(squares) show the cell-based (cluster-based) methods in the 
data, and the histograms show the equivalent distributions for 
the Gaussian noise model (see text). 



Using EJ^'^^^ it was possible to spot, in 2008, a high 
coherent noise due to the defective grounding of a barrel 
presampler HV cable and sporadic noise in a few preampli- 
fiers. These two problems were repaired prior to the 2009 
runs. The time stability of E^^^^ is regularly monitored 
using randomly triggered events by observing the mean 
and width of the E^^''''' and £'^'"" distributions. With the 
cluster-based method, the variation of all quantities was 
measured to be ±0.1 GeV over 1.5 months. This variation 
is small compared to the expected iJ™'**^ resolution (~ 5 
GeV for W ^ ev events) and can be controlled further by 
more frequent updates of the calibration constants. 

A similar analysis was performed with LI calorimeter 
triggered events, corresponding to radiative energy losses 
from cosmic muons, from the same run as used above. 
The LI calorimeter trigger (Llcalo) triggers events when 
either the sum of adjacent trigger tower transverse ener- 
gies is above 3 GeV in the EM calorimeter (EMS) or 5 
GeV when summing EM and hadronic towers [22]. The 
results are illustrated in Figure [TT] for the cell-based noise 
suppression method. Most of these events are triggered by 
energy losses in the Tile calorimeter that do not spill in 
the LAr calorimeter, which therefore mainly records noise, 
leading to a E^^^^ distribution similar to the one obtained 
with random triggers. However, in few cases, events are 
triggered by the LAr calorimeter such as the EMS trig- 
ger. The peak at 3 GeV is then shifted upwards to 6 GeV 
and the proportion of events with EJ^^^^ above 15 GeV is 
greatly enhanced. 
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Fig. 11. i?™'"" distribution with LAr calorimeter cells for 

300.000 LI calorimeter (LlCalo) triggers reconstructed with 
the cell-based method. Results for EMS trigger conditions (Sec- 
tion \3.l7S[l from the same run are superimposed on the same 
plot and the results from randomly triggered events are again 
overlaid (open symbols and histogram). 

3.3 LAr calorimeter timing 

The energy reconstruction in each cell relies on the fact 
that in the standard (five samples) physics data acquisi- 
tion mode, the third sample is located close to the signal 
maximum: this implies an alignment of the timing of all 
calorimeter cells to within a few ns. 

Several parameters determine each cell timing: the first 
contribution comes from FEB internal delays which induce 
a cell timing variation of ±2 ns within each FEB. This is 
accounted for when computing the optimal filtering coef- 
ficients. The second contribution concerns FEB to FEB 
variations due to different cable lengths to reach a given 
FEB: this relative FEB timing can vary by up to ±10 ns 
and can be corrected for by setting an adjustable delay on 
each FEB. 

The study presented here aims at predicting (using cal- 
ibration data and additional hardware inputs) and mea- 
suring (using cosmic muons and beam splash data) this 
relative FEB timing in order to derive timing alignment 
delays for each FEB. 

3.3.1 Timing prediction 

The time of the signal maximum is different in a calibra- 
tion run (tcaiib) and in a physics run (iphys)- The main 
contribution to this time is the delay Tq before the pulse 
starts to rise (the difference between the calibration and 
physics pulse widths is much smaller than this Tq delay 
variation) . This delay is driven by cable lengths which are 
different in these two configurations and additional delays 
in physics runs because of the particle time of flight, and 
the Timing, Trigger and Control (TTC) system configu- 
rations. 

In a calibration run, a signal is injected from the cali- 
bration board through the calibration cables, and is then 



read out through the physics signal cables. The value of 
the delay Tq^^^^ with respect to the signal injection can 
thus be computed for each FEB using the various cable 
lengths (icaiib, ^physics) and signal propagation speeds 

(^calib; '^physics)- 

jncalib _ -^calib ^ -^physics 
^calib ^physics 

The above prediction is compared with the measured value 
in calibration runs. The measurement corresponds to the 
time at which the calibration pulse exceeds three standard 
deviations above the noise; it is found to agree with the 
prediction to within ±2 ns, ignoring the variations within 
each FEB. 

The time of the signal maximum icaiib is obtained by 
fitting the peak of the pulse of cells in a given FEB with a 
third order polynomial. As the cable length is a function 
of the cell position along the beam axis (z, 77), the cell 
times are averaged per FEBs in a given layer (except for 
the HEC where layers are mixed inside a FEB) and a given 
rj-hin in order to align the FEBs in time. 

The time of the ionization pulse in each cell can then 
be predicted from the calibration time using the following 
formula: 

tphys = icalib — — + ^flight + ^^TTC , (8) 

l^calib 

where tcaiib was defined in the previous paragraph; ^flight 
is the time of flight of an incident particle from the in- 
teraction point to the cell, which varies from 5 ns for a 
presampler cell at 77 = 0, to 19 ns for a back cell in the 
HEC; and AtTTC is a global correction for the six par- 
titions due to the cabling of the TTC system which is 
needed to align all FEBs at the crate level. This predicted 
ionization pulse time is compared with the corresponding 
measurement in the next section. 



3.3.2 Timing measurement 

The ionization pulse time has been measured in beam 
splash and cosmic muon events. The time is reconstructed 
using optimal flltering coefficients. Since the arrival time of 
the particle is not known, one does not know in advance to 
which samples the time OFCs bi should be applied (since 
these OFCs were computed for a particular set of sam- 
ples around the pulse maximum). Therefore, an iterative 
procedure is used until the obtained At (see Eq. |4]) is less 
than 3 ns. 

The time is then corrected for two effects: first, the 
time-of-fiight difference between the beam splash or cos- 
mic muon configurations and the collision configuration, 
and second, the asynchronicity of the beam splash and cos- 
mic muon events, where arrival times vary with respect to 
the TTC clock. 

The comparison between the measured and the pre- 
dicted (Eq. |S]) ionization pulse time is shown in Figure [TH 
for the C-side (77 < 0) of each LAr sub-detector. 
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Fig. 12. Relative predicted and measured FEB times in the electromagnetic barrel (top left), electromagnetic endcap (top right), 
HEC (bottom left) and FCal (bottom right) calorimeters, for the C-side (rj < 0). The x-axis ("Slot") corresponds to a group of 
FEBs in a given layer (or a group of layers in the HEC) and rj-range. The error bars show the width of the distributions in each 
slot. 



This comparison is performed for each "sfot" corre- 
sponding to a group of FEBs in a given layer and ry-range, 
averaged over all calorimeter modules over (p. As men- 
tionned in the introduction, the relative timing of each 
group of FEBs varies by ±10 ns due to the different cor- 
responding cable lengthes. 



On the plots, the error bars correspond to the RMS 
of values for all modules in a slot: in the FCal, there is 
only one module per slot, so no error bars are shown (also 
note that slot 8 is empty in the FCal). In some regions, 
the cosmic data statistics was not sufficient to extract the 
time: the corresponding bins are thus empty. The agree- 
ment between the prediction and the two measurements 
is within ±2 ns (and at worst ±5 ns for two slots of the 
FCal). 



3.4 Signal reconstruction studies and impact on 
intrinsic global energy resolution constant term 

The main ingredient for accurate energy and time recon- 
struction of signals from LHC collisions is the prediction 
of the ionization signal shape, from which the optimal fil- 
tering coefficients used in Eq. ([3]) are computed. After re- 
calling the basics of the method used to predict the shape 
in Section [3.4.11 an estimate of the signal prediction qual- 
ity with three samples in the EM calorimeter is presented 
in Section 13.4.21 The full 32 samples shape prediction is 
used to determine the ionization electron drift time needed 
for the OFC computation in the EM calorimeter (Sec- 
tion 13.4.31) . Finally, from these two studies an estimate of 
the main contributions to the constant term in the global 
energy resolution of the EM calorimeter is given in Sec- 
tion 



Finally, a set of FEB timing alignment delays is ob- 
tained from these well understood measured relative times. 
These delays will be used at the LHC startup and updated 
once the phase between the beam and the machine clock 
will be measured and shown to be stable. The desired pre- 
cision of ±1 ns should be reached then. 



3.4.1 Prediction of the ionization pulse shape 

The standard ATLAS method for prediction of the ion- 
ization pulse shape in the EM and the HEC relies on the 
calibration system. A precisely known calibration signal 
is sent through the same path as seen by the ionization 
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pulses thus probing the actual electrical and readout prop- 
erties of each calorimeter cell. In both the EM and the 
HEC, the calibration pulse properties are parameterized 
using two variables, /stop and Tcaii, which have been mea- 
sured for all calibration boards [53] and are routinely ex- 
tracted from calibration signals [55]. 

The predicted ionization shapes are calculated from 
the calibration pulses by modeling each readout cell as a 
resonant RLC circuit, where C is the cell capacitance, L 
the inductive path of ionization signal, and R the contact 
resistance between the cell electrode and the readout line. 
The effective LC and RC have been estimated from a fre- 
quency analysis of the output calibration pulse shape |28) . 
They were also measured with a network analyzer during 
the long validation period of the three cryostats 27,'3D1 
[3T] . For the HEC, calibration pulses are transformed into 
ionization signal predictions using a semi-analytical model 
of the readout electronics, with a functional form with ze- 
ros and poles accounting for the cable and pre-amplifier 
transfer functions [35[[33] . The prediction of both the EM 
and HEC ionization pulses requires the knowledge of the 
electron drift time in liquid argon (Tdrift), which can be 
inferred from the calorimeter properties or directly mea- 
sured from data (see Section [3.4.3|) . 

To illustrate the good quality of the pulse shape pre- 
diction, radiating cosmic muons depositing few GeV in a 
cell have been used. Figure [13] shows a typical 32-sample 
pulse recorded in the barrel (top left) and the endcap (top 
right) of the EM calorimeter, as well as in the HEC (bot- 
tom left). In each case, the pulse shape prediction, scaled 
to the measured cell energy, agrees at the few percent level 
with the measured pulse. 

As already mentioned, in the FCal the calibration pulse 
is injected at the base-plane of the front-end crates, and 
therefore the response to a calibration signal differs signifi- 
cantly from the response to an ionization pulse, preventing 
the use of methods described above. Instead, seven sample 
pulse shapes recorded during the beam test campaign [9] 
110) have been averaged to obtain a normalized reference 
pulse shape for each layer. Figure [T3] (bottom right) shows 
a typical example where the agreement between the refer- 
ence pulse shape and the data is at the 4% level. 

3.4.2 Quality of signal reconstruction in the EM calorimeter 

Several PeV were deposited in the full calorimeter in LHC 
beam splash events. As an example. Figure [HI shows the 
energy deposited in the second layer of the EM calorime- 
ter. The structure in reflects the material encountered 
by the particle flux before hitting the calorimeter, such 
as the endcap toroid. In this layer, a total of 5 x 10^ five 
sample signal shapes with at least 5 GeV of deposited en- 
ergy were recorded. These events were used to estimate 
the quality of the pulse shape prediction for every cell. 
For this purpose, a Q^-estimator is defined as : 



2008 SPLASH EVENTS EM LAYER 2 E[GeV] 




Fig. 14. Total energy deposited in the LHC beam splash events 
in every cell of the EM calorimeter second layer. Empty bins 
are due to non functioning electronics. 

where the amplitude A (Eq. ([3])) is computed with a num- 
ber of samples iVsampics = 3 (because the timing was not 
yet adjusted everywhere for the beam splash events, not 
all samples can be used) , Sj is the amplitude of each sam- 
ple j, in ADC counts, gf^^^ is the normalized predicted 
ionization shape and fc is a factor quantifying the relative 
accuracy of the amplitude A. Assuming an accuracy of 
around 1%, with the 5 GeV energy cut applied one has 
'''noise < [kA)"^ . In this regime, it is possible to fit a 
function with 3 degree of freedom on the x TVdof dis- 
tribution over cells in the central region (where the 
variation is small). Therefore, A^dof = 3. A given value of 
can be interpreted as a precision on the amplitude at 
the level kQ. 

Figure [TS] shows the Q^-estimator in the second layer 
of the EM calorimeter averaged over </>, assuming k — 
1.5% corresponding to ~ 1 for 77 ^ 0. The accuracy 
is degraded by at most a factor of ~2 (i.e. ~ 4) in 
some endcap regions. This shows that these data can be 
described with a reasonable precision. 



3.4.3 Ionization electron drift time measurement in the EM 
calorimeter 

During the 2008 cosmic runs, half a million pulses with 
32 samples were recorded in the EM calorimeter from 
cells in which at least 1 GeV was reconstructed. Given 
the good accuracy of the predicted signal undershoot (see 
Figure [T3]), the drift time can be extracted from a fit to 
the measured signal [34j . 

Figure [H] shows the fitted drift time for all selected 
cells in the second layer using the standard pulse shape 
prediction method (Section [3XT]). In the EMB, the drift 
time has also been measured with a method in which the 
shape is computed using a more analytical model and 
LC and RC extracted from network analyzer measure- 
ments [SO]. The drift times extracted from the two meth- 
ods are in excellent agreement, giving confidence in the 
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Fig. 13. Typical pulse shapes, recorded during the cosmic ray campaign, for a given cell in the second layer for the barrel (top 
left) and the endcap (top right) of the EM calorimeter, as well as in the first layer of the HEC (bottom left) and in the third 
layer of the FCal (bottom right). The relative difference between data and prediction is indicated by triangles on the right scale. 




Fig. 15. Estimator (defined m the text) as a function of rj 
for 5 X 10^ pulse shapes with E > 5 CeV in the EM calorimeter 
second layer cells. is defined in Eq. ^ with k — 1.5%. 



Fig. 16. Drift time measurement in the cells of the EM cal- 
orimeter second layer with E > 1 GeV for the 2008 cosmic 
muon run. The dots correspond to drift time values averaged 
in (j). 



results: a constant value around the expected 460 ns is 
obtained, except near the electrode edges {\ri\ = 0, 0.8 and 
1.4) where the electric field is lower. The decrease of the 
drift time in the EM endcap (1.5 < \ri\ < 2.5) reflects 
the decrease of the gap size with |?7|. Similar results are 
obtained for the first and third layers of the EM calorime- 
ters. 



3.4.4 Impact on the global energy resolution constant term 
of the EM calorimeter 

When five of the production EM calorimeter modules were 
tested individually in electron beams, the global constant 
term c of the energy resolution formula was measured to 
be c ~ 0.5% in the EM barrel and 0.7% in the EM end- 
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cap [3]. The main contributors are the signal reconstruc- 
tion accuracy, the LAr gap uniformity, and the electronics 
calibration system. The first two contributions csr and 
Cgap can be investigated using results presented in Sec- 
tion 13.4.21 and 13.4.31 considering only the second layer of 
the EM calorimeter where most of the electromagnetic 
shower energy is deposited. 

From Figure [T5] . one finds that < >^ 1.4 in the 
EM barrel and 2.6 in endcap, and hence < k >— 1.8% 
and 2.4% respectively. This corresponds to residuals be- 
tween the predicted and measured pulses of 1 to 2% of the 
pulse amplitude (see Figure [T3] for illustration) , for sam- 
ples around the signal maximum. Similar residuals were 
obtained in the electron beam test analysis [55]. At this 
time, the contribution of the signal reconstruction to the 
constant term was estimated to be csr = 0.25%. Given 
the measured accuracy with beam splash events, the beam 
test result seems to be reachable with LHC collisions. 

The drift time measured in Section 13.4.31 is a function 
of the gap thickness (wgap) and the high voltage (V): 

TdHft - -f^ (10) 

where a ~ 0.3 is empirically determined from measure- 
ments [12]. In the EM barrel, the electric field is constant, 
except in transition regions, and thus the drift time unifor- 
mity directly measures the LAr gap variations. To reduce 
statistical fluctuations, the measured drift time values are 
averaged over regions of Ar] x Z\0 = 0.1 x 0.1. The distri- 
bution of the average drift time is shown in Figure [17] for 
the second layer of the EM barrel calorimeter. 



amplitude applying a factor a/{l + a) to the above re- 
sult. Therefore, the drift time uniformity leads to a dis- 
persion of response due to the barrel calorimeter gap vari- 
ations of (0.29t^;|5^)% where the systematic uncertainties 
are included. This represents an upper bound on the cor- 
responding constant term Cgap. 

For comparison, during the EM calorimeter barrel mod- 
ule construction, the LAr gap thickness was measured, 
yielding an estimate of the constant term due to gap size 
variations of Cgap = 0.16% [14] • The measurement of the 
gap size uniformity presented here takes into account fur- 
ther effects like deformations in the assembled wheels and 
possible systematic uncertainties from the in situ cosmic 
muon analysis. 

4 In situ EM calorimeter performance with 
cosmic muons 

In the previous sections, we demonstrated the good per- 
formance of the electronics operation and the good un- 
derstanding of the energy reconstruction. The cosmic ray 
events can therefore now be used to validate the Monte 
Carlo simulation that will be used for the first collisions. 

Two such analyses are presented in this section: the 
first study aims to investigate the electromagnetic bar- 
rel calorimeter uniformity using ionization signals from 
quasi-projective cosmic muons, and the second aims to re- 
construct electromagnetic showers from radiative cosmic 
muons and to compare the measured shower shapes with 
simulation. 
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Fig. 17. Distribution of the local average drift time values in 
ArjxA(j> = 0.1x0.1 bins, for the middle layer of the EM barrel. 

The drift time uniformity, estimated as the ratio of the 
RMS of this distribution to its mean value, is 1.28±0.03%. 
Using the relation between the drift time and the gap from 
Eq. [To] and the fact that the signal amplitude is propor- 
tional to the initial ionization current (/ ~ p^b<^v ~ wT" 

^ -i drift S^P 

where p is the linear density of charge), one can relate 
the relative variation of the drift time to the one of the 



4.1 Monte Carlo simulation 

The ATLAS Monte Carlo [35''36' simulates the interaction 
of particles produced during LHC collisions or from cosmic 
muons within the ATLAS sub-detectors. It is based on the 
Geant4 toolkit [37| that provides the physics lists, geome- 
try description and tracking tools. For cosmic muons, the 
material between the ground level and the ATLAS cavern 
is also simulated, i.e. the overburden and the two access 
shafts. The simulated cosmic ray spectrum corresponds 
to what was measured at sea level p8]. Air showers are 
not simulated but have a negligible effect on the analyses 
presented here. In order to save CPU time, the generated 
events are filtered before entering the full Geant4 simula- 
tion by requiring that the particles cross a specific detector 
volume (in the following analyses, typically inner detector 
volumes). 

An important use of the simulation, amongst many 
others, is to validate the selection criteria on shower-shape 
for high-level trigger and offline algorithms, as well as to 
derive the electron and photon energy calibrations. 

It is important to note that, thanks to the digitiza- 
tion step of the calorimeter simulation which emulates 
the behavior of the electronics, the standard energy re- 
construction procedure can be applied to the simulated 
events. The special procedure used for asynchronous cos- 
mic muon data, which uses an iterative determination of 
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the event time, is however not apphed to the Monte Carlo 
data. 



4.2 Uniformity of the electromagnetic barrel 
calorimeter 

4.2.1 Goals and means of the analysis 

Any non-uniformity in the response of the calorimeter has 
a direct impact on the constant term in the energy reso- 
lution (see Section [3. 4. 4p : great care was taken during the 
construction to limit all sources of non-uniformity to the 
minimum achievable, aiming for a global constant term be- 
low 0.7%. The defauh ATLAS Monte Carlo simulation em- 
ulates the effect of the constant term, but for the present 
analysis, this emulation was turned off. 

The uniformity of the calorimeter was measured for 
three barrel production modules using electrons during 
beam test campaigns |4,. Cosmic muons provide a unique 
opportunity to measure the calorimeter uniformity in situ 
over a larger number of modules, unfortunately limited 
to the barrel calorimeter due to both the topology of the 
cosmic muon events and the choice of triggers. The scope 
of this analysis is nevertheless quite different than in the 
beam test. First, muons behave very differently from elec- 
trons: in most events, they deposit only a minimum ion- 
ization energy in the liquid argon and they are much less 
sensitive to upstream material. The result can therefore 
not be easily extrapolated to the electron and photon re- 
sponse. Second, the cosmic run statistics are limited, so 
uniformity cannot be studied with cell-level granularity. 
The goal of this cosmic muon analysis is rather to quan- 
tify the agreement between data and Monte Carlo, and to 
exclude the presence of any significant non-uniformity in 
the calorimeter response. 

A previous uniformity analysis using cosmic muons [39| 
from 2006 and 2007 relied on the hadronic Tile calorim- 
eter to trigger events and to measure the muon sample 
purity. For the 2008 data discussed here, both the muon 
spectrometer and inner detector were operating and were 
used for triggering and event selection. The data sample 
consists of filtered events requiring a reconstructed track 
in the inner detector with at least one hit in the silicon 
tracker. The tracks are also selected to be reasonably pro- 
jective by requiring that their transverse (|do|) f^nd lon- 
gitudinal (l^ol) impact parameters, with respect to the 
center of the coordinate system be smaller than 300 mm. 

4.2.2 Signal reconstruction 

In the first step, a muon track is reconstructed in the 
inner detector. For that purpose, a dedicated algorithm 
looks for a single track crossing both the top and bottom 
hemispheres. This single track is then extrapolated both 
downward and upward into the calorimeter. 

Around the two track impact positions in the calorim- 
eter, a rectangle of cells (the cell road) is selected in the 
first and second layers (the signal to noise ratio for muons 



is too low in the third layer). The cells of the first layer 
have a size of Arj x Acj) — 0.003 x 0.1 and 12 x 3 such cells 
are kept. Similarly, the cells of the second layer have a size 
of Ar] X A(j) = 0.025 x 0.025, and 5x5 such cells are kept. 

To reconstruct the energy of the selected cells, the 
muon timing is obtained via an iterative procedure that is 
usually only applied to cells with an ADC signal at least 
four times the noise level. Since most muons are minimum 
ionizing particles, the muon signal is small, typically 150 
MeV is deposited in the most energetic cell in the sec- 
ond layer, only five times the noise, and many cells do not 
pass this threshold. Therefore, an alternative reconstruc- 
tion is used in this analysis: in the first pass, the iteration 
threshold is lowered to zero so that the timing is com- 
puted for most of the cells. In the second pass, the timing 
of the most energetic cell determined in the first pass is 
applied to all the other cells of the road. The cell energy is 
reconstructed at the electron energy scale and thus does 
not represent the true energy loss of the muon. Finally, 
clusters are formed in each layer to reconstruct the muon 
energy loss. The criteria used to decide on the cluster size 
are described below. 



4.2.3 Optimization of the uniformity measurement 

In order to perform the most accurate evaluation of the 
calorimeter uniformity, the measurement granularity, the 
cluster size and the selection cuts have been optimized. 
The granularity chosen is a compromise between the need 
for high statistics (large binning) and the need for high 
precision. The cluster size optimizes the signal to noise ra- 
tio while the selection cuts reduce the biases while keeping 
high statistics. 

The binning is determined by requiring a minimum of 
500 events per unit. In the rj direction, this corresponds 
to bins of 0.025 (equal to the second layer cell width) up 
to \ri\ = 0.7 and wider bins above. 

In the first layer, the muon energy loss is measured 
using a Arj x Z\(/) = 2 x 1 (in first layer cell unit) cluster, 
which contains most of the deposited energy. Adding an 
additional cell brings more noise than signal. In the second 
layer, a 1 x 3 (in second layer cell unit) cluster is used: it 
suffers less from noise than a 3 x 3 cluster, but requires 
the removal of non-projective events which leak outside 
the cluster along the rj direction. 

This projectivity cut is based on the centrality of the 
muon in the second layer cell: when the muon passes close 
to the edge of the cell, a very small non-projectivity in- 
duces a large energy leakage into the neighboring cell. 
Therefore, for each second layer cell, eight bins corre- 
sponding to the eight first layer cells located in front of it 
were defined, and in each bin a cut is applied on the beam 
impact parameter zq of the track, such that the muon is 
geometrically contained in the second layer cell. The re- 
maining statistics after this projectivity cut is 76 k events 
in the data sample and 113 k events in the Monte Carlo 
sample. The events are mainly located under the cavern 
shafts leading to a coverage of around 20% of the full elec- 
tromagnetic barrel calorimeter. 
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A comparison of the energy reconstructed in the first 
and second layers between data and Monte Carlo events 
is shown in Figure 1181 Because the muon energy loss is 
mostly ry-dependent, both distributions are shown for all 
events (top), showing a large width due to the variation 
of the energy response over 77, and for a single ij-hm (bot- 
tom). 
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Fig. 18. Energy m a 2 x 1 cluster m the first layer (histogram 
for Monte Carlo and triangles for data) and in a 1 x 3 cluster 
in the second layer (histogram for Monte Carlo and full circles 
for data) for all events (top) and a single rj-bm (bottom). 



The agreement between the data and Monte Carlo dis- 
tributions is very good, both for the shape and for the ab- 
solute energy scale which differs by only 2% in the front 
layer and 1% in the second layer. Part of the difference 
comes from the slight difference in acceptance for data 
and Monte Carlo, as well as from the difference in energy 
reconstruction. This overall energy scale difference is cor- 
rected for in the MC in the rest of the study. 



4.2.4 Calorimeter uniformity along 77 

Given the limited statistics of the projective cosmic muon 
data, the uniformity of the response in 77 cannot be esti- 
mated at the cell level. A natural choice of cell combina- 



tion is to integrate clusters in (/) since the response should 
not vary along this direction due to the symmetry of the 
calorimeter. The response along the 77 direction for cosmic 
muons depends on the variation of the amount of liquid 
argon seen by the muon. In particular, a transition occurs 
at 1 77 1 — 0.8 where the lead thickness goes from 1.53 mm 
to 1.13 mm. 

The estimation of the muon energy in each rj-hm is 
done with a fit of the cluster energy distribution using a 
Landau function convoluted with a Gaussian. The Landau 
function accounts for fiuctuations of the energy deposition 
in the ionization process and the Gaussian accounts for 
the electronic noise and possible remaining fluctuations. In 
particular, a 10% difference is observed between the width 
of the Gaussian expected from the electronic noise and the 
width of the fitted Gaussian. Mostly this bias comes from 
remaining cluster non-containment effects which are found 
to be 77-independent and thus do not produce any artificial 
non-uniformity. The most probable value (MPV) of the 
Landau distribution estimates the energy deposition. 

Distributions of data and Monte Carlo MPVs along 
the 77 direction for the first and second layers are shown 
in Figure [T^ 
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Fig. 19. Landau MPV as a function of r) in the first (top) and 
second (bottom) layers for the data (red points) and Monte 
Carlo (grey bands). 



In the first layer, the MPVs are roughly constant along 
t;, except around 7; = where some cells are physically 
missing in the detector, and around |7;| = 0.6 where the 
cell depth is varying. In the second layer, the response 
follows a typical "V-shape" corresponding to the varia- 
tion of the cell depth along 77 that rises up to |7;| = 0.6. 
Again, the agreement between the data and Monte Carlo 
is very good, showing that the contribution of systematic 
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effects due to the energy reconstruction method or the 
non-pro jectivity of the tracks is small. 

The response uniformity J7,noas is given by the RMS 
of the normalized differences between the data and Monte 
Carlo MPVs in each 77-bin : 



2008 COSMIC MUONS 



EM BARREL LAYER 1 



(11) 



- ATLAS 




; i^pVd^,^ — 

^xp. local Stat, error _ 
X 1% reference — 




r\\\ 1 ' 


1 1 ' 1 J; 2x1 clusters _ 



with: 



A^PT/, Data - MPV, 



2008 COSMIC MUONS 



EM BARREL LAYER 2 



MPF,.Data 



(12) 



where J/i^meas is averaged over 0, A^b is the number of 
bins in ry, and <C/i,mcas>=0 since the global energy scale 
difference was corrected by rescaling the MC. 

The measured uniformity should be compared to the 
expected uniformity t/cxp, which is obtained similarly to 
Eq. [TT]with ?7i,oxp given by: 



with: 



MPV, 



i,MC 



Data 



z,Data(MC) 



g(MPT/j^Data(MC)) 
-Af-PV^i,Data(MC) 



(13) 



(14) 



where cr(MPVi. Data(MC)) is the statistical uncertainty on 
the measured Landau MPV. This uncertainty is due to 
the finite statistics of the data and Monte Carlo samples 
in each bin, the Landau dispersion of the ionization, and 
the electronic noise. 

The measured uniformity t/moas should agree with the 
expected uniformity C/exp if the Monte Carlo simulation re- 
produces the data well: the key ingredients are the accep- 
tance, the muon spectrum, and the energy reconstruction 
method. A significant departure of the measured unifor- 
mity from the expected one would be a measurement of 
additional non-uniformities /7a [U\ = U^^^^ — ?7oxp)- 

The measured and expected uniformities for the two 
EM layers are shown in Figure [501 

The fluctuations of the measured energies are large: 
the RMS of the corresponding distribution is 2.4 ± 0.2% 
in the first layer and 1.7 ±0.1% in the second layer, show- 
ing that the statistical power of the analysis is limited 
given the available data and Monte Carlo statistics. The 
fluctuations mostly remain within the limits of the band 
representing the expected values. The RMS of the latter 
distribution is 2.2 % in the flrst layer and 1.6 % in the 
second layer. This demonstrates that no significant ad- 
ditional non- uniformity (Ua) is present in the data. An 
upper hmit is derived and yields Ua < 1.7% @ 95% CL in 
the first layer, and Ua < 1.1% @ 95% CL in the second 
layer. 

The calorimeter response uniformity along 77 (aver- 
aged over (f)) is thus consistent at the percent level with 
the Monte Carlo simulation and shows no significant non- 
uniformity. 
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Fig. 20. Measured (7i,meas (red points) and expected C/i,exp 
(light grey band) cosmic muon energy dispersions as function 
of rj for the first (top) and second (bottom) layers of the EM 
barrel. The dark grey band indicates a ±1% strip for reference. 



4.3 Electromagnetic shower studies 

The second analysis aims at validating the Monte Carlo 
simulation of the distribution of some key calorimeter vari- 
ables used in the ATLAS electron/photon identification. 
This is done using radiative cosmic muons that can give 
rise to electromagnetic showers in the calorimeter through 
bremsstrahlung or pair conversions. 



4.3.1 Selection of radiative muons 

To increase the probability of the presence of a muon in 
the event, it is requested that at least one track has been 
reconstructed in the inner detector barrel with \do\ < 220 
mm and pt > 5 GeV: these cuts ensure a similar accep- 
tance for data and Monte Carlo. 

A radiative energy loss is searched for in the electro- 
magnetic barrel calorimeter by requiring a cluster with 
an energy greater than 5 GeV. Since the radiation can 
occur anywhere along the muon path, the corresponding 
shower is not always fully contained in the electromag- 
netic calorimeter: this is visible in Figure [H] which shows 
the fraction of the cluster energy deposited in the first 
layer for simulated single photons from interaction ver- 
tex and for electromagnetic showers from radiating cosmic 
muons. This shows that the longitudinal shower develop- 
ment of the radiative photons is well reproduced by the 
Monte Carlo simulation, and that most of the radiating 
muons deposit very little energy in the first layer. To se- 
lect "collision-like" showers, this fraction is requested to 
be greater than 0.1. A total of 1200 candidates remain in 
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the data sample and 2161 in the Monte Carlo after this 
selection. 
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Fig. 21. Fraction of cluster energy deposited in the first layer 
of the electromagnetic barrel calorimeter for cosmic data (dots) 
and Monte Carlo (rectangles), as well as for simulated single 
photons of 5 GeV momentum from interaction vertex ( red his- 
togram). 



4.3.2 Shower shape validation 

Various shower shape distributions used for photon identi- 
fication have been compared with the Monte Carlo simula- 
tion: Figures [22] and [23] show two distributions of variables 
related to lateral shower containment in the first and sec- 
ond layers of the electromagnetic calorimeter. 
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calorimeter given by the ratio of the energy deposited in a 3x7 
cluster to a 7 x7 cluster for radiative cosmic muon data (dots) 
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Figure [22] shows the ratio of the energy deposited in 
a Ar] X A(t> = 3 x 7 (in second layer cell unit) cluster to 
that in a 7 X 7 cluster, in the second layer of the barrel 
calorimeter. In LHC collisions, this variable distinguishes 
electromagnetic showers, contained in 3 cells in 77, from 
hadronic showers, leaking outside these 3 cells. The con- 
tribution from the noise explains that the ratio can be 
above 1. 

Figure [23] shows the variable F^idc = {E±3 — E±i) / E±i 
computed as the ratio of energy within seven central cells 
in the first layer (i?±3), outside a core of three central 
cells {E±i), over energy in the three central cells : in LHC 
collisions, this variable typically separates photons, where 
little energy is deposited outside the core region, from 7r°s, 
where the two photons produced by the tt'^ deposit some 
energy outside the core region. The agreement between 
the Monte Carlo simulation and the cosmic ray data is 
very good in both the cases where the electromagnetic 
shower develops in the "collision-like" direction (in the 
bottom hemisphere) and the case where it develops in the 
backward direction (in the top hemisphere). 
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Within the statistics available from data, important 
calorimeter variables used in the electron/photon identi- 
fication in ATLAS illustrate the good agreement between 
the Monte Carlo simulation and electromagnetic showers 
from radiative cosmic events in the calorimeter. These 
results, as well as the numerous comparisons done with 
beam test data 2,3,4.5,6 , give confidence that robust 
photon and electron identification will be available for 
early data at the LHC. 



5 Conclusions and Perspectives 

The liquid argon calorimeter has been installed, connected 
and fully readout since the beginning of 2008. Since then, 
much experience has been gained in operating the system. 
Thanks to the very stable cryogenics and electronics oper- 
ation over this period, first performance studies with the 
complete LAr calorimeter coverage have been done using 
several months of cosmic muon data and with LHC beam 
splash events from September 2008. These data provided 
a check of the first level trigger energy computation and 
the timing of the electronics. In the EM calorimeter, de- 
tailed studies of the signal shape predictions allow to check 
that, within the accuracy of the analysis, there is no extra 
contribution to the dominant contributions to the intrinsic 
constant term of the energy resolution. This indicates that 
the reach of a global constant term of 0.7% is achievable. 
The non- uniformity of the EM barrel calorimeter response 
to cosmic muons is consistent at the percent level with the 
simulated response. Finally, the electromagnetic shower 
profiles are in good agreement with the simulated ones, 
thus validating the Monte Carlo description. All these re- 
sults allow for strong confidence in the readiness of the 
LAr calorimeter for the first LHC collisions. 

The ultimate LAr calorimeter performance will be as- 
sessed with collision data: this is particularly true for the 
electromagnetic and hadronic energy scale computation in 
the ATLAS environment, which is needed for many AT- 
LAS physics analyses. 
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